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Abstract: We discuss various aspects of multi-instanton configurations in generic multi-cut 
matrix models. Explicit formulae are presented in the two-cut case and, in particular, we obtain 
general formulae for multi-instanton amplitudes in the one-cut matrix model case as a degen- 
eration of the two-cut case. These formulae show that the instanton gas is ultra-dilute, due 
to the repulsion among the matrix model eigenvalues. We exemplify and test our general re- 
sults in the cubic matrix model, where multi-instanton amplitudes can be also computed with 
orthogonal polynomials. As an application, we derive general expressions for multi-instanton 
contributions in two-dimensional quantum gravity, verifying them by computing the instanton 
corrections to the string equation. The resulting amplitudes can be interpreted as regularized 
partition functions for multiple ZZ-branes, which take into full account their back-reaction on 
the target geometry. Finally, we also derive structural properties of the trans-series solution to 
the Painleve I equation. 
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1. Introduction 



The computation of instanton effects in generic matrix models is an important problem which has 
many applications. For example, matrix models describe, in the so-called double sealing limit, 
noncritical or minimal (super)string theories, and nonperturbative effects in the matrix model are 
related to the nonperturbative structure of these string theories (see [1] for an excellent review). 
Indeed, these instanton effects have been interpreted in terms of D-brane configurations (see for 
example [2, 3]) and were also instrumental for the discovery of D-branes in critical string theory 
[4]. Such connections motivated the first computations of nonperturbative effects in matrix 
models, in the pioneering works of David and of Shenker [5, 6, 7, 8], and recent progress in 
understanding the continuum side [9] has further stimulated the analysis of instanton effects also 
on the matrix model side (such as, for example, [10, 11]). 

However, most of the progress in the study of these effects in matrix models has focused on 
the double-scaling limit, i.e., when the matrix model is in the vicinity of a critical point. The 
study of nonperturbative effects for generic values of the coupling constants is also very important 
per se, and recently it has been further stimulated by the discovery that matrix models, away 
from criticality, describe in some cases topological string theory on certain non-compact Calabi- 
Yau threefolds [12, 13, 14]. With this motivation in mind, two different approaches have been 
very recently developed with the goal of computing instanton effects in generic matrix models 
off-criticality. The first approach, presented in [15], builds on previous computations in the 
double-scaling limit [6, 10, 16], extending them away from criticality. It is based on saddle-point 
techniques, and essentially computes the instanton amplitTides directly in terms of a saddle-point 
integral. The result can be expressed in terms of generating functions of correlation functions in 
the perturbative theory. Results for the one-instanton amplitude at two loops, in matrix models 
with a single cut, were obtained with this technique in [15], but going to higher order in instanton 
number and/or the number of loops seems rather complicated with this technique. The second 
approach, presented in [17], is based on the method of orthogonal polynomials, and generalizes the 
techniques first introduced in [18] in order to deal with instanton calculus. A clear advantage of 
this approach with respect to the latter is that it allows for a clean evaluation of multi-instanton 
amplitudes. However, since not all phases of matrix models can be appropriately described with 
orthogonal polynomials (the multi-cut case being a notorious example), the methods of [17] can 
not be universally applied, although they are computationally very powerful. As such, both [15] 
and [17] yield results only for matrix models with a single cut. 

In this paper we address multi-instanton configurations in multi-cut matrix models, hoping 
to go one step further with respect to [15, 17]. It turns out that, for a generic multi-cut matrix 
model background with fixed filling fractions {i.e., with a fixed number of eigenvalues in the cuts), 
a multi-instanton configuration is just any other choice of filling fractions. Formal expressions for 
the instanton amplitudes can then be obtained just by expanding the partition function around 
a generic filling fraction, making in this way contact with the results of [19, 20]. Therefore, 
the problem of obtaining multi-instanton amplitudes in the generic multi-cut matrix model is 
solved in principle by the problem of computing the partition function for a generic multi-cut 
background, and the complexity of evaluating nonperturbative effects reduces to the complexity 
of computing such a generic partition function. 

Herein, we shall consider generic multi-cut backgrounds where the eigenvalues can sit both 
at local minima as well as local maxima, which are the relevant configurations when addressing 
topological strings. Such a generic, unstable vacuum has tachyonic directions, therefore nonper- 
turbative contributions to the perturbative partition function can be exponentially suppressed 
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or exponentially enhanced. However, we shall see that it is still possible to make sense of the 
instanton expansion of the free energy by using the properties of the saddle-point 1/N expansion 
of the matrix model, and in particular of the theta-function representation discussed in [19, 20]. 

In our approach, multi-instanton amplitudes in the one-cut model can be obtained as par- 
ticular degeneration limits of the multi-cut case, by expanding around a background in which 
all but one cuts are empty. The resulting amplitude is, however, singular and it has to be ap- 
propriately regularized. We propose a natural regularization which leads to a general formula 
for multi-instanton amplitudes in one-cut matrix models. We then perform various tests of 
this formula in two examples: the cubic matrix model and its double-scaling limit, namely, 
two-dimensional quantum gravity. In both cases multi-instanton effects can be independently 
computed by using trans-series expansions of the recursion relations defining the free energy. In 
2d gravity, this recursion takes the form of the famous Painlcve I equation for the specific heat. 
The results agree with our general formula, therefore justifying our regularization. Notice that 
the appearance of divergences in multi-instanton amplitudes has been noticed before, in the par- 
ticular case of the double-scaling limit, where they correspond to singularities in the amplitudes 
for multiple, coincident ZZ-branes [21, 11]. However, it is also known that these divergences 
should be absent once the back-reaction of the ZZ-branes is appropriately taken into account. 
Our regularization procedure can then be interpreted as the correct prescription to incorporate 
this back-reaction and, in particular, our results concerning multi-instantons in 2d gravity are 
regularized amplitudes for multiple ZZ-branes in the (2, 3) model coupled to gravity. 

The organization of this paper is as follows. We begin in section 2 by reviewing both the 
general structure of multi-cut matrix models and the definition of the multi-instanton configu- 
rations we shall address. The general partition function of the matrix model is defined as a sum 
over all possible filling fractions and this leads to the tlicta-function representation of the par- 
tition function, which we relate to the multi-instanton story. In the following section 3 we then 
review the solution of multi-cut matrix models, in the large N limit, with a special emphasis on 
presenting explicit formulae for the two-cut case. In particular, we compute, in the two-cut case, 
all quantities that we earlier addressed in section 2. In this way, we are ready to shift gears in 
section 4 and consider multi-instantons in the one-cut matrix model as a particular, degenerate 
case of our previous two-cut considerations. We compute multi-instanton amplitudes up to two 
loops and discuss the nature of the instanton gas in the matrix model. We then address some 
examples, such as the cubic matrix model and 2d quantum gravity. Both cases may also be tack- 
led by making use of trans-series techniques, which provide an extra check on our results. We 
end with a concluding section, and also collect some more technical details in a few appendices. 

2. Multi-instantons and Multi— Cut Matrix Models: General Aspects 

In this section we discuss general aspects of multi-cut matrix models and multi-instantons. We 
first present a quick review of multi-cut matrix models, following [7, 19, 22, 23] (more technical 
results will be discussed in the next section). Then, we relate multi-cut matrix models to multi- 
instanton configurations in a general setting. 

2.1 Multi-Cut Matrix Models 

We recall that the one-matrix model partition function is 




(2.1) 
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Figure 1: An example of a potential leading to a three-cut solution. The N eigenvalues split into three 
sets, Ni, I — 1,2,3, and they sit around the extrcma of the potential. The support of the density of 
eigenvalues is the union of the three intervals Ai, I ~ 1, 2, 3. 



where V{x) is a potential which we shall take to be a polynomial. This can be written in diagonal 
gauge in terms of eigenvalues as 

1=1 \ i=l / 

where A (A) is the Vandermonde determinant. As usual, we are interested in studying the model 
at large but keeping the 't Hooft coupling 

t = Ng, (2.3) 

fixed. In this limit the model is described by the density of eigenvalues 

p{X) = 1^(Tt5{\-M)). (2.4) 

Let us assume that the potential V{x) has s extrema. The most general saddle-point solution 
of the model, at large A^, will be characterized by a density of eigenvalues supported on a disjoint 
union of s intervals 

s 

C=[jAi, (2.5) 

7=1 

where Aj = [x2/_i, X27] are the s cuts and xi < X2 < ■ ■ ■ < X2s- If s > 1 we call this saddle-point 
a multi~cut solution of the Hermitian matrix model. A potential leading to a three-cut solution 
is depicted in Fig. 1. This multi-cut saddle-point can be described in terms of an integration 
over eigenvalues as follows. In the s~cut configuration described above the A^ eigenvalues split 
into s sets of Nj eigenvalues, / = 1, . . . , s. Let us denote each of these s sets by 

{a1'^K=i,...,7V„ I = 1,...,s. (2.6) 

The eigenvalues in the I-th set sit in the interval Ai around the /-th extremum. Along this 
interval, the effective potential 

yeff(A) = V{\) - t J dX'p{X') log |A - A'l (2.7) 
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is constant. It is possible to choose s integration contours Cj in the complex plane, I = 1, . . . , s, 
going off to infinity along directions where the integrand decays exponentially, and in such a way 
that each of them passes through exactly one of the s critical points (see for example [24]). The 
resulting matrix integral is then convergent and can be written as 



N 



The overall combinatorial factor in this expression, as compared to the one in (2.1), is due to the 
fact that there are 

iV' 

ways to choose the s sets of Ni eigenvalues. Of course, when the integrand is written out in 
detail, it splits into s sets of eigenvalues which interact among them through the Vandermonde 
determinant (see for example [25]). In this paper we will assume that the filling fractions, 

ei = ^ = J^dXpiX), 7 = 1,2,. (2.10) 

are fixed, and they can be regarded as parameters, or moduli, of the model. One may also use 
as moduli the partial 't Hooft couplings tj = tej = ggNj, which can be written as 



with ^2^1=1 — Here y{z) is the spectral curve of the matrix model which, in the multi-cut 
case, describes a hyperelliptic geometry, and relates to the effective potential as V^^{z) = y{z). 
We shaU be more concrete about it in the next section. Notice that, in general, the saddle-points 
we consider will have unstable directions. This is in fact the generic situation one finds in the 
applications of matrix models to topological string theory [12, 13, 14]. Indeed, as we shall learn 
in this paper, precisely the fact that we consider general saddle-points will make it possible to 
extract multi-instanton amplitudes from multi-cut matrix models. 

The free energy of the multi-cut matrix model has a genus expansion of the form 

oo 

F = log Z = Y,Fg{ti)gl^-\ (2.12) 

9=0 

In particular, the planar free energy -Fo(^/) can be computed from the spectral curve by using 
the special geometry relation 



dFo{t) 
dti 



I dAy(A), / = l,...,s, (2.13) 

Jb{ 



where is a path which goes from the endpoint of the Ai cycle to the cut-off point A. This 

point may then be taken to infinity after one removes any divergent pieces from the integral. The 
higher Fg{ti) can also be computed explicitly, and there is a systematic, recursive procedure to 
obtain them, put forward in [26, 27]. 
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2.2 Multi-Instantons in Matrix Models 

If one regards (2.8) as the matrix integral in a specific topological sector, characterized by the 
fillings Ni, . . . , Ns, it is then natural to consider the general partition function as a sum over all 
possible arrangements of eigenvalues across the several cuts [5, 7, 19, 20] 

^(Ci,---,C.)= E Cf •••C^^(iVi,...,iV,). (2.14) 

Ni + -+Ns=N 

The coefficients may be regarded as ^-parameters which lead to different ^-vacua, and we 
refer the reader to [7] for such an interpretation. The sum (2.14) may also be regarded as a 
matrix integral where the N eigenvalues are integrated along the contour 

s 

C = Y,CkCk\ (2.15) 

k=l 

in this case the ^-parameters yield the relative weight of the different contours C^, see [7, 20]. 

A very important point is the fact that the different sectors appearing in (2.14) can be re- 
garded as instanton sectors of the multi-cut matrix model. Indeed, let Z{Ni, . . . ,Ns) be the 
partition function in one sector with filling fractions {Nj}, and Z{N[,--- ,N'^) be the parti- 
tion function corresponding to a different choice of filling fractions {Nj}. In this case one can 
immediately write 

This expression means that, if we pick a set of filling fractions (A'^i, . . . ,Ns) as our reference 
background, all the other sectors will not be seen in Qg perturbation theory. It is then only 
natural to regard them as different instanton sectors of the matrix model. 

Furthermore, notice that, depending on the value of the real part of the exponent appearing 
in (2.16), the sectors with filling fractions {N[, ■ ■ ■ , N'g) will be either exponentially suppressed 
or exponentially enhanced, with respect to the reference configuration (iVi,--- ,Ns). Let us 
consider for example a cubic potential where A'^i eigenvalues sit at the minimum of the potential 
and A''2 eigenvalues sit at the maximum. The sector with fillings (A^i — 1,A^2 + 1), which may 
be regarded as the one-instanton sector, will be more unstable than our reference configuration, 
and typically it will be exponentially suppressed. However, the anti-instanton sector with fillings 
(A'^i + 1,A^2 — 1) is more stable and it will be exponentially enhanced. Therefore, in multi-cut 
matrix models (and in contrast to most field theories) instantons and anti-instantons have actions 
with opposite signs. As we shall see in the following, this leads to some subtleties. 

Let us then consider a reference configuration characterized by arbitrary filling fractions 
{Nj}, I = 1, . . . , s, and regard any other set of filling fractions {N'^} as instanton configurations 
around our reference background. In this case, the determination of the instanton expansion of 
the partition function amounts to a Taylor expansion around the configuration with tj = Qs^j. 
For simplicity we will write down formulae for the two-cut case, and refer to [20] for the formal, 
general expressions. We start from a reference configuration {Ni,N2) and write 



Ni Ni 

C^z{Ni-e,N2 + e) = z(''\Ni,N2) C^^^^^ (2-i7) 

e=-N2 e=-N2 
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where we have denoted zW(A^i, A?2) = Z{Ni,N2) and 

^ ZiN^-i,N, + i) 

Zi0){N„N,) ■ ^'-'^^ 

Let us now assume that Nj ioi I = 1, 2. We recall that ti = ggNj and, as such, in terms of 
the 't Hooft parameters one may equivalently write 

vffl _ Z{ti-£gs,t2+£gs) 

" Zio)it,,t2) • ^'-'^^ 

In the two-cut case it turns out to be rather convenient to use the variables 

t = ti + t2 and s = ^(ti-t2). (2.20) 

In this case, one immediately finds, by simply expanding around = (or, equivalently, around 
s = 0) 

= g4 exp |l - g,[i^sF^ + ^ d^Fo) + 0(5.^)}. (2.21) 

In this equation we have set 

A{ti) = dsFoiti) and g = exp(a2^o)- (2.22) 

The first thing one recognizes is that A corresponds to the action of an instanton obtained by 
eigenvalue tunneling from the first cut to the second cut. Notice that, due to (2.13), we can write 
A as an integral over the spectral curve, obtaining the familiar result 



A= dxy{x). (2.23) 



The expression (2.21) gives a general formula for the first terms in the Qs expansion of the 
^-th instanton contribution to the partition function, in a general two-cut background. Notice 
however that, since £ can be either positive or negative, the expansion in (2.17) will be a Laurent 
expansion in both ^ = exp{—A/gs) and In this case, the free energy will also be expressed 
as a Laurent scries in ^ and but each coefficient in this scries will be given by an infinite 
sum of terms. This is the subtlety we were previously referring to, and it comes from the fact 
that we are expanding around saddle-points which may have unstable directions. 

There is a rather natural way to re-sum these terms by using theta-functions, as first in- 
troduced in [19, 20]. At large Ni and N2 one can extend the sum in (2.17) from —00 to +00, 
obtaining in this way 

Z = Z^^\t^,t2) C'e^^ 1-5. + - 93^0)+ 0(5,2) . (2.24) 

If one now exchanges the sum over £ with the expansion in gg, we may immediately write Z in 
terms of the Jacobi theta-function 

+00 2 

Mq\z) = 1^ (2.25) 
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where z = ^ 



Z = Z(0)(ti,t2) {MqIz) - 9s{zd,MQ\z) d.Fi + i (zd^fM^lz) d^Fo) + Oig^,)\. (2.26) 



6 

As we shall explicitly show in the next section, this theta-function is well defined due to the fact 
that 

^ (2.27) 
27ri ^ ' 

is indeed the modulus of an elliptic curve, hence 1^1 < 1. One may now use Jacobi's triple identity 

oo oo oo 

M<l\z) = 11(1 - 11(1 + zg-V2) + ^-ig-V2) (2.28) 



n=l n=l n=l 



to write 



where 



oo / 1 \£ e , —i 

log UqIz) = log m + E ^ ^r^^> (2-29) 

e=l * q2 - q 2 



<t>{q) = ll{l-q^). (2.30) 

n=l 

The last two terms may also be written in terms of the quantum dilogarithm 



e=l * q2 - q 2 

This reorganization makes it possible to express the total free energy F = log Z in terms of an 
infinite series which formally has the structure of an instanton/anti-instanton expansion, as 

F = ih, t2) + log 0(g) + C' e-'^/^^ (l + Oigs)) . (2.32) 

e^o ^(«2 - g 2 ) ^ ^ 

It is of course possible to write the Qs corrections in a similar way. Notice that log 0(g) gives a 
contribution to -Fi coming from instanton/anti-instanton interactions in the partition function. 

In this paper we shall mostly explore another avenue. We will expand the multi-cut partition 
function (2.14) around the most stable configuration, which generically is a one-cut configuration. 
In this way we will obtain general formulae for multi-instanton contributions in the one-cut 
matrix model by using multi-cut matrix models, generalizing the one-instanton formulae we 
first obtained in [15]. Before doing that, in the next section we shall give details on the solution 
of multi-cut matrix models which will then be useful in our subsequent computations. 



3. Solving Multi-Cut Matrix Models 

In this section we will review the explicit solution of the multi-cut matrix model, focusing on the 
two— cut case. In particular, we will explicitly compute the several relevant quantities appearing 
in the multi-instanton formulae we have obtained in the previous section. An introduction to 
many of the results collected in the following may be found in, e.g., [19, 22, 23]. 
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3.1 General Results 

In order to solve for the density of eigenvalues (2.4) we introduce, as in the one-cut case, the 
resolvent 



^ ' fc=0 

which has a standard genus expansion uj{z) = Yl^=o9s^ ^gi^) 

{z) = JdX^. (3.2) 



The normalization of the eigenvalue density 



J dAp(A) = 1 (3.3) 



implies that 

Mz) - ^ (3.4) 
as z ^ +00. Notice that the genus-zero resolvent determines the eigenvalue density as 

p{z) = -r^ {uJo{z + ie) - uJo{z - ie)) . (3.5) 

ZTTl 

One may compute loo{z) by making use of the large A'^ saddle-point equations of motion of the 
matrix model, 

Looiz + ie) + cooiz - ie) = -V'{z) = 2 P / dA (3.6) 

t Jc z — X 

For a generic multi-cut solution, the large N resolvent is given by 



„„(,) = i./^^X>),M4, (3.7) 

^ ^ 2t Jc 2Tri z - w ]l a{w) ' ^ ' 

where we have introduced the notation 

2s 

a{x) = l[{x-Xk). (3.8) 

k=l 

In (3.7), C denotes a closed contour enclosing the union of intervals appearing in the standard 
multi-cut case, as we have discussed in the previous section. An equivalent way to describe the 
large N solution is via the spectral curve y{z), which is given by 

y{z) = V'{z)-2tu;o{z) = M{z)^/^, (3.9) 

where 

(3.10) 

is the moment function. The moments of M(z), which are defined by [22] 

M^'^=<f^ K>) , (3.11) 
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can be easily calculated to be 



and will play an important role in the following. We shall also denote 

Mi = mI^^ = M{xi). (3.13) 

In order to fully determine the large N solution one still needs to specify the endpoints of 

the s cuts, {xk}k=i 2s- The large z asymptotics of the genus-zero resolvent immediately yield 

s + 1 conditions for these 2s unknowns. They are 

"^^'^"^ - 2a„., (3.14) 



C 27ri 



w] 



for n = 0, 1, . . . , s. In order to fully solve the problem, one still requires s — 1 extra conditions. 
These are obtained by using (2.10) and treating these filling fractions (or, equivalcntly, the partial 
't Hooft couplings (2.11)) as parameters, or moduli, of the model. Since X]/=i ~ (2-10) 
one indeed only finds s — 1 conditions, as initially required. 

Later on we shall need explicit expressions for the derivatives of the endpoints of the cuts, 
Xj, with respect to the 't Hooft moduli, tj. These are obtained as solutions to a linear system 
which we will now write down. Taking derivatives with respect to tj in (3.14), one immediately 
obtains 

Y^MiX^-^ = A6ks, k = 0,...,s. (3.15) 

i=l ^ 

For fixed j, this gives s + 1 conditions for 2s quantities. The remaining s — 1 conditions can be 
obtained precisely from the 't Hooft moduli, 

tk = ^ n d\M{X)./\^\, (3.16) 

by taking derivatives with respect to tj. A similar calculation is done in [22], in its Appendix A. 
By deforming the contour away from infinity, we can write the moment function as 

M(A,^-/ (3.17) 

where Cx is a small contour around w = X. Therefore, by taking derivatives with respect to tj in 
(3.16), we obtain 

2s 



; P dA </ ^Il^f^ L_]^/M. (318) 

47r /cue. 27ri (a; - A) Va; - X-XiJ^^' ^ 

where the derivative acted on cr{X) and a{io). Since 



A — (J 



LO — Xi X — Xi {UJ — Xi){X — Xi)^ 



(3.19) 
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we finally obtain 

r ^p±(_^._J_)^^.:mM,. (3.20) 

JcuCx 27ri [uj - X) \u - Xi A - Xj/ -^(7(0;) A - 

Notice that the integrand no longer has a pole at a; = A, hence only the integral around C 
contributes. We then end up with the following s equations 



2s 



j=l ^'■J •JX2k-l 



If we introduce the integrals [22] 

Ki,k= / dA^P^, (3.22) 



X — X 



we may then write 

2s 



^f2MiKi,k-^ = 6jk, k = l,...,s. (3.23) 
^ i=i ^ 

These equations, together with (3.15), fully determine the derivatives dxi/dtj. Observe that 
we have written down 2s + 1 equations in this procedure, but, as it turns out, only 2s of these 
equations are actually independent, uniquely determining our 2s unknowns. Other derivatives 
with respect to the 't Hooft parameters can also be expressed in terms of the derivatives of the 
branch points. For example, from 

it is simple to find 

dtj 2 ' dtj 2^ xi-xk dtj' ^ ' ' 

3.2 The Two Cut Matrix Model 

In the two-cut case one can write very concrete results for the different quantities involved in the 
above solution, in terms of elliptic functions [22, 19, 25]. We now review some of these results. 
The support of the density of eigenvalues is given in this case by 

^1 U^2 = [a;i,X2] U [x3,a;4]. (3.26) 

We define, as in [19], 

ic = r = ' Kik), k- = (^i-^2)(^3-x4) 

\/{xi-X^){x2-Xa) {Xi-X2,){x2-Xi) 

where K{k) is the complete elliptic integral of the first kind. Similarly, 
K' = r = ' K{k% k'^ = l-k^= (3.28) 
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Wc refer the reader to the appendix for a hst of definitions and conventions concerning elliptic 
functions. One further defines the modular parameter r of the elliptic curve = (j{x) by 

It is now easy to show that the second derivative appearing in the one-loop coefficient of the 
£-instanton sector, (2.22), is given by 

5 = 27rir. (3.30) 

To do this, we adapt an argument from [19]. Prom the special geometry relation (2.13) it follows 
that 



rX3 

j dx{V'{x)-2tuQ{x)), (3.31) 



hence 



(3.32) 



Using the asymptotic behavior of ojq{x) in (3.4) it is easy to prove that 

d{tujQ{x)) _ C{s,t) 
ds 

where C{s,t) does not depend on x. Its value can be fixed by using that 



(3.33) 



ti= (f Prtiooix). (3.34) 
JAi 27ri 

Indeed, since 
it follows that 

C{s,t) = -^. (3.36) 

Now plugging this result inside (3.32), we immediately find (3.30). 

In the two-cut case, the equations for dxi/dtj which were obtained in the previous section, 
can also be explicitly written in terms of elliptic functions. The required integrals Ki^k were 
evaluated for alH = 1, . . . , 4, and for = 2, in Appendix B of [25], and we can use the expressions 
therein to solve the resulting system of equations. One may also use a procedure identical to the 
one in [25] to obtain the Ki^k integrals when k = 1, with the same final result for the dxi/dtj 
derivatives. In particular, one finds the resulting simple answer 

dx^^j7r_ 1 

ds MiK:Uj^,{xi-xjy ^ • ' 

Using the above formulae, it is easy to obtain explicit expressions for the other derivatives 
of the free energy appearing in (2.21), in particular the two-loop contribution to the £— instanton 
sector. The cubic derivative dgFo can be obtained directly from (3.30) as 
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The first factor in tlie riglit liand side may be computed witli tfie use of the identities 

dp- 2P^^^ 

while the second factor is more involved and can be obtained from (3.37). The result, for a 
generic two-cut matrix model, is 

dk'^ _ 27r 1 1 

57 ~ Mi---M^ K{k) (Xi - X^fl\x2 - X4)3/2 Wi^^iXi - Xj) ' 

■ (M1M2M3 (Xi - X2fixi - ,X3)2(X2 - X3f + 

+M1M2M4 (Xi - X2)^(X1 - X4)^(X2 - X4)^ + M1M3M4 (xi - X3)2(xi - X4)^(x3 - X4)^ + 
+M2M3M4 (X2 - X3)^ (X2 - X4)^ (X3 - X4)^) . (3.40) 

Putting everything together we finally obtain 

1 d^Fp _ TT^ i/(xi -X3)(X2 -X4) 1 

6 ds^ 6M1 • • • M4 K^{k) (Xi - X2)(X2 - X3)(xi - X4)(X3 - X4) lli^jiXi - Xj) 

■ (M1M2M3 (Xi - X2)^(X1 - X3)^(X2 - X3)^ + 

+M1M2M4 (Xi - X2)^(xi - X4)^(X2 - X4)^ + M1M3M4 (xi - X3)^(xi - X4)^(x3 - X4)^ + 
+M2M3M4 (X2 - X3)^(X2 - X4)^(X3 - X4)^) • (3.41) 

Another quantity which can be computed in the two-cut case in terms of elliptic functions 
is Fi(ti,i2) [22, 25]. It is given by 

4 

^^'^ = - ^ E 1^^* - ^ l^^(^) -J^H - ^j)' + I - ^3)' + I ln(x2 - X4f. (3.42) 

1=1 i<j 

This is, in fact, the last ingredient required in order to compute the complete two-loop coefficient 
in the £-instanton sector expansion, where we are still missing the expression for dgFi in the two- 
cut case. As it turns out, the resulting expression is rather long, and not particularly illuminating 
and, as such, we present it in the appendix. 

4. Multi— Instantons in the One— Cut Matrix Model 

We shall now shift gears and consider multi-instantons in the one-cut matrix model as a par- 
ticular, degenerate case of our above two-cut considerations; in this way generalizing the one- 
instanton results we have previously obtained in [15]. 

4.1 General Formulae 

As we mentioned at the end of section 2, one of our goals in this paper is to focus on the 
expansion of the multi-cut partition function around its most stable cut, suppressing in this way 
the anti-instanton configurations. This is the natural choice, for example, if we are interested in 
computing a convergent matrix integral in terms of a perturbative series in gg plus exponentially 



-13- 



Figure 2: The total partition function (4.2) is a sum over all possible instanton sectors. The ^-th 
instanton sector is obtained by tunneling £ eigenvalues from the most stable, one-cut configuration, to 
the other (unstable) saddle-point. 

small corrections, as in [17]. We will assume that we are in a generic situation and the potential 
has a unique global minimum (like in Fig. 2). If we call Ai the cut surrounding this minimum, 
the most stable configuration has filling fractions 

(iV,0,---,0) (4.1) 

and it is a one-cut solution of the model. Any other configuration will be exponentially suppressed 
with respect to this one, and a convergent multi-cut matrix integral can be computed at large 
by considering the partition function Z{N, 0, • • • ,0) and then adding exponentially suppressed 
contributions from the other configurations. In this way we obtain a multi-cut configuration 
as a sum of multi-instantons in the one-cut matrix model. In particular, this produces general 
formulae for the multi-instantons amplitudes of the one-cut solution. 

We will focus for simplicity on stable one-cut configurations which arise as limiting cases 
of two-cut matrix models. This means that we will only consider one type of instanton, arising 
from eigenvalue tunneling from the filled cut to the empty cut. The generalization to models 
with more than two cuts should be completely straightforward^ . 

The total partition function is obtained by summing over all possible numbers i of eigenvalues 
tunneling from the filled cut to the empty one, as 

N 

z = Y,C'z{N-e,e). (4.2) 

We are interested in the 't Hooft limit in which iV — > oo with t = QgN fixed, and i <^ N. From 
the point of view of the spectral curve, the partition function Z{N—i, i) corresponds to a two-cut 
curve in which one cut is filled with N — i eigenvalues, with N ^ 1, and the other cut is nearly 
but not completely pinched, since it contains i <^ N eigenvalues (see Fig. 3). In the 't Hooft 
limit the finite sum in (4.2) becomes an infinite one, and one can write 

Z = Z(0)(t)|l + f;c'^^'^|, (4.3) 
e=i ^ 

^The situation we analyze in here is similar to the "birth of a cut" transition studied in [28]. For previous work 
on phase transitions of this type in the Hermitian matrix model see, e.g., [29, 30, 5]. 
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Figure 3: The partition function in the i-th instanton sector corresponds to a spectral curve in which 
one of the cuts is filled with N — £ eigenvalues, with N ^ 1, while the other cut is nearly pinched, since 
it contains £ ^ N eigenvalues. 



where Z^^\t) = Z{N,0) is the total partition function of the one-cut matrix model, i.e., the 
model in which the eigenvalues sit at the minimum, and 

(,) _ z{t-eg„egs) 

Just like we did before, we can try to evaluate Z^^^ by expanding the numerator around gg = 0. 
If one tries this there is a subtlety, however, since the free energies Fg(ti,t2) are not analytic at 
t2 = 0. Geometrically, this corresponds to the fact that we are expanding around a configuration 
in which the second cut of the spectral curve is completely pinched. To understand the origin of 
this nonanalyticity, write 

Fg{ti,t2)=F^{t2)+Fg{ti,t2), (4.5) 

where F^{t) are the genus g free energies of the gauged Gaussian matrix model with 't Hooft 
parameter t, i.e., 

F,^{t) = ^t^{logt-l), Ff(t) = -^logt, ••• (4.6) 

It is simple to see that Fg{ti,t2) is indeed analytic at ^2 = 0, so that the lack of analyticity of the 
matrix model free energy at ^2 = is due to the Gaussian contribution F^{t2). This issue has 
been clarified in [31] in a slightly different context: the Gaussian part of the matrix model free 
energy comes from the measure, and it is not analytic when the 't Hooft parameter vanishes. The 
"regularized" Fg{ti,t2) comes from re-summing the perturbation theory double-line diagrams 
with genus g and it is an analytic function. 

Physically, the reason for the appearance of this singularity is that in this problem t2 = igs 
and i is small as compared to A^. Therefore, it is not appropriate to treat the integration over 
the i eigenvalues from the point of view of the large expansion. Instead, they should be 
integrated exactly. This argument also suggests that, in order to regularize the computation, we 
should subtract F^{t2) from the total free energy and at the same time multiply Z^^^ by the 
exact partition function Z^ for the gauged Gaussian matrix model with i eigenvalues, 

^F=(|)^G^(^ + 1)' (4-7) 
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where G2{i + 1) is the Barnes function 



-1 



(4.9) 



G2(,i + l) = Y[i\. (4.8) 

i=0 

The appropriate expression for the partition function around the ^-instanton configuration is 
then 

= exp \Y,9's'-' [Pgit - ksAds) - Fg{t)) 

^9>0 

We can now expand the exponent appearing in (4.9) around gs = 0, since it is analytic, and we 
obtain, up to two loops, 

Z^'^ = (^^2(^ + 1) C'g^ exp ^^l-g,{i^sF^{t) + J^^Jo{t))+0{9',)y (4.10) 

In this equation 

A{t) = dsFo and q = exp(d^Fo). (4.11) 



Notice that all derivatives appearing in both (4.10) and (4.11) are evaluated at ti = t and t2 = 0. 
We can now compute the free energy as 

oo 

F = logZ = ^/FW (4.12) 

i=o 

where we set 

z = Ce"^/^^ (4.13) 

Since we are expanding around the most stable configuration Z now has a Taylor series in z, 
and F will be well-defined as a formal series in positive powers of z. This series is appropriately 
written in terms of "connected" contributions 

f +00 oo 

e=i ^ e=i 



where 

z(2 = Y^tDll ^ = + . (4.15) 

s>l * fciH \-ks=i k=l 

We then deduce that the -f-instanton contribution to the free energy of a one-cut matrix model 
is given by 



The explicit expression of F^^^ at two loops is 

1/2 



/ f(^) = (^e Zi^\. (4.16) 



i{2TT)2 [ \ 6 J J 



yielding formulae for the two-loops ^-instanton amplitude in an arbitrary one-cut matrix model. 
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4.2 The Nature of the Instanton Gas 

An interesting point to notice is that our result for Z^^\ (4.10), is not typical of a dilute instanton 
gas. Indeed, in a dilute instanton gas one has that, at leading order in the coupling constant, 



(4.18) 



If this was the case here Z^^^ would scale as gi^"^, but it is clear from (4.10) that this is not the 

case: Z^^^ scales as gl ^"^ . The reason why this happens is due to the presence of the Vandermonde 
determinant, and it can be easily understood from a simple scaling argument. The ^-instanton 
integral is roughly of the form 



« / J]dAiA2(A)exp 



1 



-2^KffME(^^-^o)' 



(4.19) 



1/2 

If we change variables to Ui = (Aj — XQ)/gs it immediately follows that the measure for the (. 

£/2 

eigcnvahics leads to the factor gj typical of a dilute instanton gas. However, in the same way 
the Vandermonde determinant leads to an extra power of g^s^ The instanton gas in a matrix 
model should be rather regarded as an ultra-dilute instanton gas since, at weak coupling <C 1, 
the partition function for £ instantons is even more suppressed than in the usual instanton gas. 
Physically, the ultra-diluteness is of course due to the eigenvalue repulsion. 

In particular, we disagree with the analysis of multi-instantons in the one-cut matrix model 
proposed in section 2 of [10]. In that paper, the Vandermonde interaction between the instantons 
is set to one and the resulting integral factorized. This cannot be done without jeopardizing the 
very scaling of Z^^^ with ^r^.. We shall later on verify this scaling, as well as other details of our 
main formulae (4.10) and (4.17), against other techniques for computing multi-instanton effects. 

4.3 Explicit Calculations 

Our main multi-instanton formulae (4.10) and (4.17) involve varioTis derivatives of two-cut free 
energies, in the limit in which one of the cuts shrinks down to zero size, i.e., in the i2 limit. 
We will now derive explicit expressions for these derivatives and, in particular, we shall check 
that the expression obtained for F^^^ agrees with the result we previously obtained in [15]. 

The limit ^2 — ^ corresponds geometrically to a degeneration of the curve in which the A2 
cycle shrinks to zero size, i.e., ^ X/^ = xq. The spectral curve of the two-cut problem becomes 



y{x) M{x){x - xq)^{x - xi){x - X2). (4.20) 

Notice that here M{x) is the moment function of the original two-cut problem. The moment 
function of the one-cut problem, obtained by degeneration, is then given by^ 

Mi(x) = M{x){x - xq). (4.21) 

We can now compute the quantities appearing in the ^-instanton free energies (4.17). First 
of all, the instanton action is given by 



A{t) = ]im (dsFoit, t2) - dt,F^{t2)) • (4.22) 



^Recall that, in the one-cut problem, it is precisely the requirement M\(xo) = that allows for the nontrivial 
saddle-point xo characterizing the one-instanton sector [15]. 
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Both quantities appearing here are regular at t2 = and, from (2.13), one immediately finds 



^ rxo 

A= dx Mi{x)y/ {x - xi){x - X2), (4.23) 

J X2 



1X2 

which is the instanton action in the one-cut case. 
We next compute 

\dlm,Q) = \\\ra (dlFo{t,h)-dlF^{t2)) = lim -\\ogt^. (4.24) 

In the limit 2:4 x-j, the elliptic modulus appearing in (3.27) vanishes and both r and logt2 
diverge. In order to calculate this limit, as well as similar ones, we will set 

e = X4 - X3. (4.25) 

As e — >^ the leading behavior of the elliptic modulus is given by 

7rir~log— , (4.26) 

i.e., it diverges as loge. According to our general argument above, this divergence should be 
removed by subtracting (logt2)/2. In an appendix we compute the expansion of t2 in powers of 
e. In particular, we find that at leading order t2 ^ e^, so indeed the cancellation takes place. 
Putting together (4.26) and (C.6) we find 

-1/2 k"^ Xi- X2 1 

q^l'^ = hm = = , (4.27) 

^^ M{xo)[{Xi- Xq){X2- Xq)\2 

Since 

M{xg) = M[{xq) (4.28) 

(4.27) is in complete agreement with the formula obtained in [15] for the one-loop contribution 
to (the factor l/v^ appearing in the formula of [15] is already included in (4.17)). 

We now proceed with the calculation of the quantities appearing at two loops. We shall first 
calculate dgF^, which is defined by 

\dlFo{t, 0) = lim Q53Fo(t, t2) - (4.29) 

and should be regular. Indeed, it is easy to see that (3.41) has a pole of order two in e = X4 — X3. 
Making use of (C.6) one can explicitly check that both the divergent pieces in e"^ and e"^ are 
removed by subtracting 1/(6*2), and we finally get 



6 * ' i/(xo - x\){xq - X2) I {xi - X2)^ixo - X2Y M{X2) (xx- X2Y{xq- xxY M{xx) 

4M'(xo)2 3M"(xo) 17(,-ri+,r2-2.ro).l/'(.r„) , 



+ 



+ 



M(xo)3 2M(xo)2 2(xo-xi)(xo-X2)M(xo)2 

llx\ + 118x1x2 + 77x^ - 272x0 (xi + X2 - xq) 
8(xo -xi)2(xo -X2)2M(xo) 



(4.30) 
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A similar computation can be done in the case of Fi. Again, from the exphcit expression 
(B.4) it is easy to see that dgFi is singular as e ^ 0, but 



= lim(a.F,(M2) + ^) (4.31) 



turns out to be regular, as expected, and can be computed to be 

gPun^ 1 f M'jxo)^ _ M"{xo) _ {xi+X2-2xo)M'{xo) 

''^'^ ^/{xo-xi){xo-X2) l6M(xo)3 8M(xo)2 24(xo - xi)(xo - ^2) M(xo)2^ 

^ 19{Xi + X2) — 22xiX2 — 16Xo (Xi + X2 - Xq) (Xq - X2)(X2 - 3X1 + 2xo) 



96(xo - Xi)2(x2 - xo)2M(xo) 3(xi - X2)2(xo - Xi)2M(xi) 

(xo - xi)(xi - 3x2 + 2xo) (xo - X2) Af'(.xi) 



+ 



3(xi - X2)2(xo - X2)2M(X2) 4(xi - X2)(xo - Xi) M(xi)2 

(xo-xi)M^(x2) 
4(xi - X2)(xo - X2) M(X2) 



(4.32) 

This result, together with (4.30), gives a completely explicit result for the amplitudes (4.10) and 
(4.17) in terms of geometrical data of the spectral curve for the the one— cut matrix model. 
From the above expressions one easily checks that 

^d^,Fo{t,0)+dsFi{t,0) (4.33) 

agrees with the two-loop result in equation (3.63) of [15], after taking into account that .xi, X2 were 
denoted by a, b in there, and that the moment functions are related by (4.21) (in an appendix, 
we include some further relations between the moment functions of the one and two-cut models, 
in order to make the comparison more explicit). This tests our expression (4.17) for at two 
loops against our calculation of the one-instanton amplitude in [15], with complete agreement. 

4.4 Multi— Instantons in the Cubic Matrix Model 

Our results (4.10) and (4.17) provide general expressions for multi-instanton amplitudes in one- 
cut matrix models. We shall now exemplify and test these formulae, within the multi-instanton 
set-up, in the cubic matrix model. To perform the test, wc have to compute multi-instanton 
amplitudes independently and, in order to accomplish this, we will use the trans-series formalism 
based on orthogonal polynomials which was recently developed in [17]. 

Let us briefly summarize the formalism of [17], to which we refer the reader for further details 
and examples. We first recall that the orthogonal polynomials Pn(A) for the potential V{X) are 
defined by 

dA -YW 

— e 3s Pn{X)PmW = hnSnm, n > 0, (4.34) 



/ 



27r 



where pn are normalized by requiring that ~ A" + • • • . It is well known (see for example [1]) 
that the partition function (2.1) of the matrix model can be expressed as 



N-l N 



Zn = U h, = llr^-\ (4.35) 



i=Q i=l 
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The coefficients 

rn = (4.36) 
satisfy recursion relations depending on the shape of the potential. They also obviously satisfy 

Tn = (4.37) 

In the limit N ^ oo, ^ becomes a continuous variable that we will denote by z, and r„ is 
promoted to a function, R{z,gs), for which we have the multi-instanton trans-series ansatz 

oo 

R{z,gs) = Y.C'R^'\z,gs)e-'^^'^/3% (4.38) 

1=0 

where the £-instanton coefficient can be written as 

R^'\z,gs) = Rf{z) (^l + f;5."i?W(z)^ (4.39) 

and C is a nonperturbative ambiguity. The coefficients Rn\z,gs) can be computed from the 
so-called pre-string equation, a difference eqiiation which may be derived as the continuum limit 
of the recursion relations satisfied by the coefficients r„. For any polynomial potential the pre- 
string equation can be written down explicitly [1]. Furthermore, once (4.38) has been found, we 
can then extract the free energy F{z,gs) = logZ from the continuum limit of expression (4.37), 
namely 

F{z + gs, gs) + F{z - gs, gs) - 2F{z, 5,) = log R{z, 5,). (4.40) 
The total free energy F{z,gs) also has an instanton trans-series expansion of the form 

Fiz,gs) = ^e-'^^'y^^F^'\z,gs), (4.41) 

where 

{z, gs) = F^'^ {z)(l + f; g: (z) j , i>l. (4.42) 

The different coefficients Fn\z) can be easily obtained from A{z) and I^\z,gs), see [17] for 
some explicit expressions. Notice that in here the nonperturbative ambiguity C is an integration 
constant and it cannot be fixed within this method. In particular, it might differ from the 
nonperturbative ambiguity C appearing in the partition function (2.14) by a multiplicative factor 
independent of z, the 't Hooft parameter (see [17] for a detailed discussion). 

We now use the above trans-series formalism to compute the amplitudes F^^^ for a cubic 
matrix model with potential 

V{z) = -z+^-. (4.43) 

A similar computation has been performed in [17] for the quartic matrix model. The recursion 
relation for the coefficients r„ reads [1] 

Tn (a/1 - rn - rn+1 + ^Jl-rn-Vn-l^ = gsU. (4.44) 
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The continuum limit of equation (4.44) above is 

R{t, gs) (^V'^-R{t,gs)-R{t + 9s,gs) + V^- R(t,gs) - R{t- gs,gs)) = t. (4.45) 

At lowest order in £ and gs we find 



2R^^\t)^l-2R^^\t)=t. (4.46) 

It turns out to be convenient to express everything in terms of a new variable, R'^\t) = r. 
Power series expanding equation (4.45) , and solving it recursively, one can compute in this way 
the coefficients Rn\r{t)). The first few read 

Rio)(t) = (9r-5) (0). . ^ 3(l62r3 + 1017r^-1316r + 385) 

2 ^' 32(1 -3r)4' 4 ^ 2048(3r-l)9 



v'SF^' 192(1 -3r) 



5 • 

2 



With these functions in hand, we can now compute the free energies using (4.40), up to an 
overall constant for each F^^\ In turn, this can be compared to our prediction (4.17). We 
find 

4(3r-l)4 I 192r(l-3r)2 

^(2-) r -8 - 228r + 429r2 oA 

F^^^ = -^\^+9s -^ + 0{gl)\, (4.48) 

32(3r-l)2 I 96r(l-3r)2 J ^ 

F (3. = -— |i + + Obi: 

192(3r-l)x I 64r(l-3r)2 

These results precisely have the structure predicted by (4.17). This can be explicitly checked by 
computing the quantities dgFi, d^Fo and q = e'^«^" from our analytical expressions (4.27), (4.30) 
and (4.32), specialized to the cubic model. We find 

r 

Q = r> 

32(1 - 3r)5 

1^,-:- -183r2 + 84r + 8 

^9tFo = (4.49) 

6 96r(l-3r)2 ^ 

-57r2 + 60r - 8 



dsFi 



192r(l - 3ry. 



If we now plug these in (4.17) we immediately recover (4.48), up to an r-independent normal- 
ization of C^. This is a rather strong check of our multi-instanton formulae. 

4.5 Multi— Instantons in 2d Quantum Gravity 

Unfortunately, the results we have just obtained for the cubic matrix model test the structure 
(4.17) only for low instanton number £. On the other hand, it is well-known (see, e.g., [1]) that 
the cubic matrix model has a double-scaling limit which describes pure 2d gravity. We shall now 
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derive analytic expressions for F^^'^ in this limit, up to two loops, and for all f > 1, by using the 
string equation, and we will verify that they are in perfect agreement with (4.17). 
The critical point of the cubic matrix model occurs at 

.0=^, r. = l. (4.50) 

The double-scaling limit is defined as 

5 5 5 1 

5s — 0, z ^ Zc, K* = {zc — z)iSsgg fixed. (4.51) 
In this limit, the specific heat of 2d gravity 

u{k) = {3R{k) - 1) g7^ (4.52) 

satisfies the Painleve I equation 

u{Kf - ^u"{k) = K, (4.53) 

as can be deduced from the difference equation (4.45). The multi-instanton contributions can 
thus be obtained by finding a trans-series solution 

u{z) = J2C'^^'\^) (4-54) 
to the Painleve I equation. One obtains a set of recursion equations for the u^^\z) of the form 

+ =0, e>l. (4.55) 

fe=0 

Changing variables to 



z-i = VSx, (4.56) 



each u^^\x) has the structure 

u(^) (x) = x"^^ e-^^/^ (x), (4.57) 

where 

A = 

5 



A = ^ and e(^)(x) = l + ^«f (4.58) 



fe=i 

In particular, for ^ = 0, 1, we have 

e(°)(x) = 1 - + 0{x^), e(i)(x) = 1 - ^ + Oix^). (4.59) 

In the trans-series solution (4.54) parameterizes the nonperturbative ambiguity. In prin- 
ciple the value of u^q^ is not fixed by the recursion and it could be absorbed in (. However, if we 
choose C to be identical to the nonperturbative ambiguity appearing in the matrix model, then 
the value of u'^^ can be fixed, as we shall do in the following. 
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From e^^\x) it is possible to compute the multi-instanton corrections to the double-scaled 

<z) = -^ds(-^)> (4.60) 



free energy, i^^^ (x) , which are defined by 



and 

oo 

i^ds(^) = EC'e-'^^'^i?(^)- (4-61) 



As a function of x they have the structure 

oo 



k=l 



The exponent Xe and coefficients /^^^ can be related to the data appearing in u^^^ . We have, for 



example, 

A£ = + f^^' = -"-^ /}^^ = - I + -. (4.63) 

We can now use the equations (4.55) to obtain recursive relations for the various coefficients 
appearing in the ansatz (4.57). The recursion for is very simple, 

2 

a£ = - + ak + ae-k, k = l,...,£-l. (4.64) 

This requires to be linear in £, and we find 

£2 £ 

In order to find F^^'ix) up to two loops, we have to solve the recursion relations for both Uq ^ 
and u^^^ . Starting with Uq^ , we find the recursion 



- X: 4^ut'\ £ > 2. (4.66) 



" 2(^2 

^ ' k=l 

Using the well-known identity 

i-i 

Y^ki£-k) = -£{£^-l) (4.67) 

k=l 

one can immediately verify that 



W _ 12^ 







4^3 



5 ^ ^ 



solves (4.66). The solution (4.68) was already obtained in [32] in their study of nonperturbative 
effects in 2d gravity, and in [33] in their study of the trans-series solution to the Painleve I 
equation. The second nontrivial coefficient in the trans-series solution, u\ , satisfies the recursion 



-P-r^pl^^-fEM*-*)"?'). «>2. (4.69) 

k=l 
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It can be easily checked that 
solves (4.66). Here one uses that 

Y,k\i-k) = y'{e-i). (4.71) 
fc=i 

Notice that (4.70) does not specialize to for £ = 1. Plugging these results in (4.63) we find 

(1) 

/o " 



4 ) 

4-35 



We can then obtain, directly from the string equation, 

(4.73) 

^rw = tit^ (/«)' {i - « (g - . + o(x^)}, £ > 2. 

On the other hand, the double-scaling limit of our multi-instanton formulae (4.17) leads to 

F^X-) = C {l - (^1,1 + ^^0,3) X + 0{x^)\, 

e-i i/2 ^^•'^^^ 
(^) = ^'^^V'' { 1 - ^ (<^i,i + ^'^0,3 + gds) X + 0(x2) } , £>2. 

In these formulae, q^g, </'o,3 and (/>i^i are the double-scaling limits of q, dgFo and dgFi, respectively. 
Their values can be obtained explicitly from (4.49) and read 

1 47 17 

9ds = -^, ct>o,3 = - and 01,1 = ^. (4.75) 

Since u^^^ = (see (4.59)), we see that the explicit solution to the recursion equations, (4.73), 
agrees with the prediction (4.74) arising from our multi-cut analysis. Agreement of both formulae 
further requires 



V27r/ 

!S th 

ambiguities are the same, as we pointed out above 



This is the normalization condition that fixes the value of u^^ so that the nonperturbative 
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It is easy to solve the string equation to higher orders and we can use the results obtained 
in this way to give further tests of our general formulae. Since the £-instanton amplitude (4.10) 
gives contributions to the free energy at i{£ — l)/2 loops, we will focus on the nonperturbative 
partition function 



Z°P = l + ^z^4;^=exp 



where 



z = Ce"^/^. 

The general formula (4.10) leads, in the double-scaling limit, to 



(4.77) 



(4.78) 



(27r) 



G2{^ + 1) gf/' \l-x(l </.i,i + - <^o,3) + C(a 



(4.79) 



We shall now present a slightly different reorganization of .Z'^p, inspired by the analysis of 
[33]. Instead of considering the partition function Z^^^ at fixed instanton number and expanded 
in X, we shall consider in the terms with a fixed power of x but for all instanton numbers. 
In other words, we write 



fe>0 



where 



and 



(4.80) 



(4.81) 



(4.82) 



We are now going to extract a structural result for Zk{S,) from the general formula (4.79). This 
formula, based on the connection to the multi-cut matrix model, says that for each i the starting 
power of X is £^/2. But the total power of x in the expression above, for each fixed i, is k + £/2. 
Therefore, (4.79) says that, at fixed k, zj-^e is different from zero only for the finite number of ^'s 
such that 

£{£ - 1) 



< A;. 



(4.83) 



Let us denote by £m{k) the maximum number satisfying this bound. We have, for example, 

C(0) = 1, £m{l) = 2, £m{2) = 2, (4.84) 



and so on. We then conclude that 



.(fe) 



(4.85) 
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is a polynomial of degree iraik)- We have computed Zi^{^) from the recursion (4.55) at high k 
and found that indeed they are polynomials in ^ of the right degree. We list them up to = 6: 

^(2^ + 111) 



ZsiO = - 

zm = 

z^iO = - 

ZeiO 



192 

$, (1048C + 19299) 
24576 ' 
C (160^2 + 11705160C + 114670521) 
70778880 ' 
^ (552320^2 + 12466492352^ + 79686828333) 
18119393280 ' 
C (331932480^2 + 3646348240864^ + 17179325749341) 
1159641169920 ' 
C (102400^=^ + 15883648256640^2 + 105786318127916160^ + 395885797054644519) 



6679533138739200 

(4.86) 

The results for Zfc(^) are in full agreement with the expected structure (4.79). For example, if 
k is such that (4.83) is saturated for imik), it follows from (4.79) that the highest power of ^ in 
Zjt(^) should be given by 

qiG2{em{k) + l). (4.87) 

This is indeed the case in all the examples we have tested, and it is furthermore a nontrivial test 
of our regularization (4.9). It says that the string equation for pure 2d gravity "knows" about 
the Barnes function, and therefore about the Gaussian matrix model partition function which 
we have used to regularize the full partition function. 

There are two interesting applications of our above results for the multi-instanton amplitudes 
in 2d gravity. First, the formulae (4.74) for F^^ should yield the partition function of £ ZZ-branes 
in 2d gravity {i.e., in the (2, 3) minimal model coupled to gravity) after including their full back- 
reaction on the target geometry. As explained in [21], these amplitudes are divergent in the 
continuum theory, but one should be able to regularize them by taking into account the fact that 
the £ ZZ-branes smooth out the pinched cycle of the Riemann surface describing the model [34] . 
In section 4.1 we performed the computation of the amplitudes starting from the theory in which 
the pinched cycle is completely smoothed out, i.e., the two-cut solution. When we expanded 
naively around the background with a pinched cycle, we found that the expansion was singular. 
However, we resolved the singularity in the matrix model computation by performing an exact 
Gaussian integration of the £ eigenvalues. This led to (4.9) and eventually to (4.74), which, as 
we have seen, agrees with the explicit trans-series solution of Painleve I. In the language of [21] 
this means that our matrix model computation appropriately incorporates the back-reaction of 
the ZZ-branes, and it is in contrast with the computation of [11] where the interactions between 
identical ZZ-branes have not been taken into account. 

The second application of our result is more mathematical. The reorganization of the in- 
stanton expansion in (4.80) can, of course, also be made at the level of the free energy and the 
specific heat. Indeed, if we write 

oo 

n = x-2/5^x'=nfc(0, (4.88) 

k=0 
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it can be seen by plugging this ansatz in Painleve I that the Uk{^) satisfy differential equations 
and are given by rational functions [33] . It then follows that in the expansion of the free energy 



the fkiO ^6 also rational. But, after exponentiation, the partition function (4.80) is written in 
terms of polynomials Zk{(,), as we have argued. This is rather surprising from the point of view 
of the original Painleve I equation and seems to be a new structural result for the trans-series 
solution to Painleve I. 

5. Conclusions and Outlook 

In this paper we have studied multi-instanton configurations in multi-cut matrix models. As we 
have seen, producing general, formal expressions for their amplitudes is rather straightforward, 
although on a general background there are subtleties related to the existence of tachyonic direc- 
tions. The focus of our paper has thus been the analysis of such expressions in the one-cut limit, 
generalizing the one— instanton formulae from our earlier work [15]. This is also a limit in which 
the generic expression needs a regularization, and we have proposed a natural one which fits all 
available data on multi-instantons, in particular the multi-instanton amplitudes in 2d gravity. 

There are various avenues for future research. Our study has focused for simplicity on two- 
cut models and their one-cut limit, since this has made possible detailed tests of our approach. 
General formulae for more than two cuts can similarly be obtained in a rather straightforward 
way (closely related formulae in the general case have already been written down in [20]). It is 
also straightforward to generalize the results of section 4 and to consider multi-instantons in a 
one-cut model which is obtained as the limit of a multi— cut model with more than two cuts. In 
this way one could incorporate different types of instantons, corresponding to tunneling from the 
filled cut to the several different empty cuts. Another venue of research deals with extending our 
formulae for the (2, 3) minimal model to general minimal string theories. This would produce 
amplitudes for the interaction of an arbitrary number of ZZ-brancs, which fully consider their 
back-reaction in the target geometry, greatly enhancing our understanding of multiple brane 
interactions, albeit in the noncritical context. Naturally, this also raises the (much harder) 
question of how to derive our results directly in the continuum Liouville theory. 

What remains one of the most interesting open issues for us would be to understand and 
test the connection between the multi-instanton formulae in the multi-cut case and the large- 
order behavior of the perturbative amplitudes for a fixed multi-cut background, following our 
general analysis in section 2. A significant step towards efficiently computing the topological 
1/N expansion of multi-cut models was very recently given in [35, 36], building on [37, 38], and 
opening way to a large-order analysis, and this is a subject we hope to return in the future. This 
would generalize the one-cut study performed in [15] and would also have direct applications in 
topological string theory on toric Calabi-Yau manifolds. 
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A. Elliptic Functions: Definitions, Integrals and Properties 

We follow the conventions in [39] . The complete elliptic integral of the first kind is defined as 



K{m) = / . (A.l) 

Jo ^/(l-^2)(l-m^2) 



V^(l-t2)(l-mt2) 

In here m = A;^ is the parameter of the complete elliptic integral and k is the elliptic modulus. 
One also defines the complementary modulus as /c'^ = 1 — k'^. The complete elliptic integral of 
the second kind is defined as 

Eim) = At (A.2) 

Both elliptic integrals, K{ra) and E{m), have a branch cut discontinuity in the complex m plane 
running from 1 to oo. The complete elliptic integral of the third kind is defined as 

U(n\m)= / , . (A.3) 

In here we take < A;^ < 1 and n is the elliptic characteristic. 

The complete elliptic integrals of first and second kind satisfy the Legendre relation 

E{k)K{k') + E{k')K{k) - K{k)K{k') = |. (A.4) 

One sometimes writes K{k') = K'{k) and E[k') = E'{k), and denotes these functions as com- 
plementary elliptic integrals. The complete elliptic integrals of first and second kind also relate 
to each other via derivation, e.g., 



dK _ E{k) - k'^K{k) 
'dk ~ kk^ ' 

dE _ E{k) - K{k) 
'dk ~ k ■ 

Useful properties satisfied by the complete elliptic integral of the third kind include 

U2 



(A.5) 
(A.6) 



, s nl-k^^fk^-n\ A;^ 

n(nm) + - n m\=-^ K{m), (A.7) 

1 — nk'^ — n \l — n J k^ — n 

n{n\m) + u( — \m]=K{m) + ^J- (A.8) 



n\ J ^ ' 2]! (l-n)(n- A;2)" 
Let us next consider the elliptic geometry given by the curve y^{z) = M'^{z) cr{z), with 

4 

Cr{z) =Y\_{Z - Xi) = {Z - Xl){z - X2){Z - Xz){z - X4) (A. 9) 

i=l 

and xi < X2 < xs < X4. A non-constant single-valued function f{z), with z G C, is called an 
elliptic function if it has two periods 2uJi and 2lo2, i.e., 

f{z + 2muji+2nuj2) = f{z), (A.IO) 
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with m, n G Z. The half-period ratio of the ehiptic geometry, r, is defined as r = and the 
eUiptic nome is defined as q = e^'^'^. The eUiptic nome can be related to the complete elliptic 
integral of the first kind as 

making it clear that q{m) has a branch cut discontinuity in the complex m plane running from 
1 to cx). The inverse relation expresses the complete elliptic integral of the first kind in terms of 
a Jacobi theta function as (here '&3{q) = i^siqlz = 0)) 

K{k) = ^^l{q). (A.12) 

Useful integrals which are used in the main text, and which can be obtained from the tables 
in [39], are 

K{k), (A.13) 



XI ^/(7{z) V(^4-^2)(^3-^ 

^4 dz 2i 



1x3 ^A(^ V^(X4 - X2){X3 - Xi] 

along the A-cycles; as well as 



K{k), (A.14) 



J X2 



\/cr{z) \/{x4: - X2){X3 - Xi) 

along the i?-cycle. In the formulae above the elliptic modulus is 



K{k'), (A.15) 



e = (^1 - ^2X0:3 -^.4) 

[Xi - X3){X2 - Xi) 



B. Expression for dsFi in the Two— Cut Matrix Model 

A final ingredient required in the main body of the paper in order to explicitly write all multi- 
instanton formulae in the two-cut matrix model, in particular in order to compute the two-loops 
coefficient, is the derivative of genus one free-energy, Fi(ti,t2)- This was computed in [25], with 
the result 

4 

^1 = E 1^ - ^ l^^(^) - ]^ E 1^ - ^^•)' + ^ 1^ - ^3)' + ^ In (^2 - x^f , (B.l) 

i=l i<j 

where we now need to evaluate 

dFi 1 1 dMi E{k) - k'^K{k) dk^ 1 ^-^ 1 fdxi dx 



ds 24^Mi~57 4k'^k'^K{k) 'ds 6 ^ Xi - xj \ ds ~ ~ds , 

1 / dxi \ 1 / dx2 dx4 

4 (xi — X3) \ ds ds J 4 (x2 — X4) \ ds ds 



(B.2) 
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The only novelty is which may be computed with some of the formulae we presented in the 
main body of the paper. The result is 

1 ^ 1 dMi 7ry(xi-X3)(x2-X4) 1 fM^ , ^ 

24 g 1^1)7 = WW) n.<,fe-x,) i w '"^ - """"^ - - 

M" M" 

(^1 - X3){xi - X4)(xs - X4) + {Xi - X2)(xi - X4){x2 - X4) - 

-J^ixi-X2){xi-X3)(x2-X3)j. (B.3) 

Assembling all different pieces together, it finally follows 



dFi TT ^ixi-X3){x2-X4) ( M'l 

- " Y72 ~ ^3){X2 - X4){X3 - X4)- 



ds 8K{k) U^<j{x^-XJ) 

M" M" 

(^1 - Xz){xi - X4){xz - X4) + (Xi - X2)(xi - X4)(X2 - X4) - 

M'' \ , T:{E{k)-k'^K{k)) (xi-X3)3/2(x2-X4)3/2 

(Xl - X2)(X1 - X3)(X2 - X3} + 



Mr ^ ' ' 2K2(A;) W.^,{x,-x,f 

TT (^2 - X3)^(X2 - X4)^(X3 - X4)^ + TT i^'^ ~ Xzf{xi - X4)^(x3 - X4)^ + 
1V1\ IVI2 

+ (Xi - X2)^(xi - X4)^(X2 - X4)^ + TT i^'^ ~ X2)^(xi - X3)^(x2 - Xs)^ ) + 
M3 M4 J 



^7r\/(;i;i - .r:-i)(.i'2 - .ri) / 3.r j - XyX2 - x^x^ - X2X1 - -ktix-^ + 2x2X->, + 2,r3X4 



^K{k) \ (Xi - X2)2(X1 - X3)2(X1 - X4)2Mi 

Co — X1X9 — xiX'i — X9X^ — 4x9X4 + 2xiX4 + 2x^X4 ^ 



(xi 


- X2)^(X2 - X3^ 


|2(X2 - 


X4)2M2 


3x| — X2X3 - 


- X2X4 — X3X4 — 


4xiX3 


+ 2xiX2 + 2xiX4 


(xi 


- X3)2(X2 - X3^ 


|2(X3 - 


X4)2M3 


8X4 — X1X3 - 


- X1X4 — X3X4 - 


- 4X2X4 


+ 2xiX2 + 2x2X3 


(xi 


- X4)2(X2 - X4^ 


)2(X3 - 


X4)2M4 



+ 



C. Series Expansion of the Vanishing Period t2 



(B.4) 



In the two-cut matrix model, the vanishing period, i.e., the period of the differential y(x)dx 
which vanishes as X4 ^ X3, is given by the following integral: 

t2 = — / dxM(x)x^^. (C.l) 

271" Jxz 

We want to find a systematic series expansion of this quantity in powers of 

e = X4 — X3. (C2) 
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This is easily done by writing 

(X - Xi){x - X2) = (X3 - Xi){X4 - X2) (l + ^—^) (l + ^ ~ ""M . (C.3) 

If we now introduce the variable 



u = (C.4) 



the integral defining t2 becomes 

t2 =— M{X3) [(X3 - Xi)(x3 - X2 + e)] 2 



JO 



, Mfx'ij + eu) / u \ 2 u — 1 

dn ■ / 1 . - 111,- 



M(X3) V 2:3 -xi/ V e + X3-X2/ ^ ^ ^ 



(C.5) 



This expression admits a systematic expansion in e. In order to go to two loops in the computation 
of multi-instanton effects, we need to compute t2 to order e^. One finds, after setting X3 = xq, 



N /7 T7 e/^ 1 1 2M'(xo) 

t2 = — M{xo)V{xo-Xi){xo-X2) 1 + T I + + — 



4 Vxo - Xi Xo - X2 M{xo) ^ 
5eV ixi-X2f ^ 4ixi+X2-2xo)M'{xo) ^M'^x^Y 



128 \(xo - xi)2(xo - X2)2 (xo - xi)(xo - X2)M(xo) M(xo) 

D. Relating Moment Functions of One and Two— Cut Models 

In this appendix we wish to relate two-cut moments to one-cut moments, following the line 
which we outlined in the main body of the paper. In the one-cut limit, one makes use of 
M\{z) = M{z){z — Xo) and may write, for z ^ xq 

M(.) = (D.1) 

Xq- Z 

f . Mi{z)\ , , 

M'{z) = M[{z) + , p.2 

Xq- Z \ Xq- Zj 

M"{z) = — (m^{z) + (m[{z) + , (D.3) 

Xq - Z \ Xq - Z \ Xq - Z ) ) 

alongside with M^'^^(xq) = m[^~^^\xo), n > 0, among several other similar expressions. 
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